Spin-flop transitions and spin- wave gaps in La2Cu04 
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We study the spin- wave spectrum and the spin- flop transitions in La2Cu04 in a uniform mag- 
netic field at zero temperature. Using the non-linear cr-model, we show that a field applied along 
the orthorhombic b direction leads to a two-step rotation of the staggered magnetization, first in 
the be and then in the ac plane, until the order parameter is completely aligned along the c axis. In 
contrast, for a perpendicular magnetic field, we find a conventional spin-fiop transition induced by 
the competition between the field and the interlayer coupling. A comparison with recent measure- 
ments of the field-dependence of the in-plane spin-wave gap shows a beautiful agreement between 
theory and experiments. 
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I. INTRODUCTION 



Since the discovery of high-temperature superconduc- 
tivity^ in materials derived from the La2Cu04 (LCO) 
parent compound, much progress has been made towards 
a better understanding of the magnetic properties of this 
Mott insulator. From a theoretical point of view, the puz- 
zle behind the physics of the cuprates lies in the interplay 
between the charge carriers and the spin background. A 
sound knowledge of the magnetic properties of the par- 
ent compound is therefore indispensable for a successful 
description of the doped materials. The relatively sim- 
ple structure of LCO facilitates this task considerably by 
allowing one to perform a great number of experiments. 

At low-energy, the charge degrees of freedom in 
LCO are frozen and the remaining copper spins are 
well described by an antiferromagnetic spin-1/2 Heisen- 
berg model. The field-theoretical treatment of this 
model with isotropic superexchange interactions was in- 
troduced by Chakravarty et alw^^ and proved to be 
very successful in describing the paramagnetic phase of 
LCO. However, it was noted early on that in the low- 
temperature orthorhombic phase, anisotropies due to the 
Dzyaloshinskii-Moriya (DM) and the direct-exchange in- 
teraction cannot be neglected^'^ and greatly influence the 
magnetic behavior in this phase. Quite recently, Chovan 
and Papanicolaou^ as well as Silva Neto et alii general- 
ized the non-linear cr-model (NLSM) approach to systems 
with anisotropic superexchange. Employing this formal- 
ism, we present a study of the evolution of the magnetic 
properties of LCO in the presence of a uniform magnetic 
field, which allows us to compare our findings with re- 
cent measurements of the spin- wave gaps determined by 
Raman spectroscopji^ and shed some light on the contro- 
versy concerning the nature of the spin- flop transitions. 
Although the mechanism of the spin-flop transition ob- 
served for a perpendicular (to the copper-oxide planes) 
magnetic field is well understood^'', and basically due 
to the competition between the antiferromagnetic inter- 
layer coupling and the magnetic field favoring ferromag- 
netic alignment, the precise value of the critical field is 



not yet univocally determinedSiii. We hope that the 
present work encourages new neutron scattering experi- 
ments that can clarify this issue. Even more controversial 
is the case of a uniform field applied along the easy axis, 
i.e., along the direction of the staggered magnetization 
observed at zero field. There is no doubt that in this 
case the order parameter is rotated towards the c axis, 
but Thio et aL~^ suggested a two-step reorientation (sup- 
ported by our calculations) while Ono and co-workersi^ 
argued in favor of a single continuos rotation. 

The rest of the paper is organized as follows. In Sec.lTTI 
we introduce the non- linear cr-model describing LCO. We 
first examine this model in the absence of a magnetic 
field and discuss the strength of the interlayer coupling 
and the renormalization of the gyromagnetic ratio due 
to the spin-orbit coupling. Sec. IIIII then contains our 
main results of the evolution of the spin-wave spectrum 
in the presence of uniform magnetic fields applied along 
the three different crystal axes. These results allow us to 
proceed to a careful comparison between our findings and 
recent experimental observations of one-magnon Raman 
spectra^ in Sec. IIVI Finally, we present our conclusions 
in Sec. El 



II. MODEL 

The NLSM is a very convenient tool to describe the 
low-energy dynamics of the Heisenberg model, relevant 
for LCO. In this approach, the staggered component of 
the copper spins in a single layer of LCO is represented 
by a continuous vector field n {r) of unit length n{v)^ — 
1. To avoid confusion, we denote vectors acting in the 
three-dimensional (3D) spin space by arrows and vectors 
acting in the 2D coordinate space by the usual bold font. 
We consider a system of L layers and indicate a given 
layer by a subindex I, i.e. (r). Throughout this paper, 
we adopt the orthorhombic coordinate system, with unit 
vectors Cq,, a = a,b,c. Coordinate and spin space are 
linked through the pinning of the commensurate Neel 
magnetization to the orthorhombic b axis. It is therefore 
convenient to use the same coordinate system in both 
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cases with real space unit vectors = . According to 
Refs. 7 8, the Lagrangian reads 



Hi (fii X + (ni X 



l} , (1) 



- 2B ■ [D X ni] - Yni ■ m+i 



with the "hard" constrainti4 n(r)^ = 1. Here x± ~ 
« 4.5 • 10^^ {meV)~ is the magnetic susceptibihty 
of the Heisenberg model and c — 1.17V2J « 230 meV 
the spin-wave velocity (u sing the exchange coupling J « 
140 meV given in Ref. Ila) . Note that due to quan- 
tum fluctuations in the Heisenberg model, X-L is re- 
duced by 50% compared to its classical value. The 
anisotropies are due to the DM interaction, with a DM 
vector directed along the orthorhombic a direction of 
length D « 2.5 meV, and the XY-term, which leads 
to ~ 5 meV. These values follow from neutron 

scatteringi^. The magnetic field B is given in units of 
gUB = g 5.79 • 10~^ meV/T, g being the gyromagnetic 
ratio and fiB the Bohr magneton. The coupling between 
spins in adjacent copper-oxide layers is described by the 
last term in Eq. Q ■ The coupling constant 7^ can be de- 
termined via the spin-wave dispersion (see Sec. Ill A|) or 
the critical field of the spin-flop transition (see Sec. lHl C|l . 



A. Interlayer coupling constant 

Let us first study the Lagrangian ^ in the absence of 
a magnetic field and determine the spin-wave dispersion 
in order to obtain an estimate of the strength of the inter- 
layer coupling constant 7^. In the static limit, the stag- 
gered magnetization is oriented along the orthorhombic b 
direction (due to the anisotropies) and we can therefore 
consider excitations of the form 



ni (r) 



(r) 



"i,2 (r) 



For small deviations from this ground state, the time 
derivative is given by 

ni (r) = 

\ni,2 (r)^ 

Introducing the Fourier transformation 



{ql+k-l 



(2) 



and substituting these expressions in the Lagrangian Q , 
we find 



C 



EE 

i=l,2 q 



i2nf 



-n 



with = and M| = F^. Since the prefactor x±/2 
and any terms independent of Hi (r) are not important 
for the purpose of the present work, we omit them from 
now on. The equations of motion read 



c-'k' + +Y {I + cos q) 



resulting in the following spin-wave spectrum 

ujj c^k'^ + +-f'^{l + cos q) , 
Lu'^ = c^k^ + Fc + 7^ (1 + cos q) , 

with the usual time dependence 



"k,9j (t) = 



(3) 



(4) 



Both, the in- and out-of-planc modes are gapped, due to 
the presence of the anisotropies Mj. Comparing Eas. 
with the spin- wave dispersion given in Eq. (2) of Ref. llSl 
we find 



7^ = 2a±c^ w 5 (meVy 



(5) 



using a± = 5 • 10^^. Although a comparison with the 
experimentally determined spin-wave spectrum is theo- 
retically possible, we would like to point out that the 
above value of a± is actually based on estimates of the 
spin-flop transition presented in Ref. and, as far as 
we understand, does not follow from neutron scattering 
experiments. 



B. Gyromagnetic ratio 

The assumption that the copper spins in LCO are ac- 
curately described by the usual electron gyromagnetic 
ratio g = 2 is not quite correct. Let us therefore consider 
the renormalization of g due to the spin-orbit interac- 
tion. The electronic configuration of the Cu in La2Cu04 
is with one active hole per Copper atom. In the case 
of non-zero orbital angular momentum I, the spin-orbit 
coupling and the interaction with the magnetic field read 

Hir^t^Xl-S-fiB {l + 2S^B , 

with A « —0.09 eV, see e.g., Ref. ^3 Considering the 
terms involving the orbital angular momentum as pertur- 
bations, we derive an effective Hamiltonian of the form 

Hctf — —ga f^B S^B" , 

with a generalized gyromagnetic tensor g which is diag- 
onal in the axes a, 6, and c. The ground state and the 
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TABLE I: Lowest-lying energy levels and wave functions of 
LaaCuOi. For simplicity, we neglect the small orthorhombic 
distortion. Yi^m are spherical harmonics. 



Orbital Wavefunction Energy (eV) 

M^2_y2 1^2,2 + ^2,-2) 

M^y 15^2,2-^2,-2) 1.35" 

3d3,2_,2 |y2,o) 1.5' 

3d,, 1^2,1 + ^2,-1) 1.70" 

idy. \Y2,i - Y2,-i) 1.70" 



In order to keep the formalism as light as possible, we 
chose to restrict the staggered magnetization to be uni- 
form within a given copper-oxide layer, i.e., ni (r) — fti. 
This approach has two advantages: Firstly, we do not 
have to carry the in-plane momenta through the deriva- 
tions, which greatly improves the readability of the pre- 
sentation and secondly, it is straightforward to remove 
this constraint, by restoring the momenta in the final 
answer, by replacing wj' ^ w| — c^k^. 

A. Magnetic field along a 

For a magnetic field applied along the a axis, i.e., B — 
Bca, the Lagrangian (Q) reads 



first three excited states of La2Cu04 are summarized in 
Tab.Q] For simplicity, we neglect the small orthorhombic 
distortion here. Second order perturbation theory yields 
corrections to the energy given by 



6E„ 



{0,a\l ■ B\n,a) {n,a\l ■ S \0,a') 
Eo — En 



where the factor 2 gives the correct multiplicity of the 
considered processes. With the knowledge of the non- 
zero matrix elements of the orbital angular momentum 
operator 

{3dxy\ l'^ \3dx2_y2'^ — 2 , 

{3dyz \ r \ 3dx2-y2) — 1 , 

{3dxz \ I'" \ Mx2_y2) = i , 
it is straightforward to calculate the gyromagnetic ratio 



9 2A 

q I =2 

^ 1.7 



511 =2- 



8A 
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2.1 , 
« 2.5 



(6) 



We use these values for further calculations. 



III. FIELD DEPENDENCE OF THE SPIN- WAVE 
SPECTRUM 

Depending on its direction, an applied magnetic field 
can lead to a reorientation of the staggered magnetiza- 
tion. Similar to the approach followed in Sec. Ill Al in the 
absence of a magnetic field, we proceed as follows. We 
first find the static ground state orientation of the or- 
der parameter and then consider excitations, described 
as small deviations from the ground state configuration. 
This allows us to determine the spin-wave spectrum in 
the presence of uniform magnetic fields. The excitations, 
denoted rii and 712 are perpendicular to the static ground 
state orientation and correspond to the in- and out-of- 
plane modes at zero magnetic field. 



L 



1=1 



-B' 



{n\) +K)' 



Clearly, a field along the a direction does not change 
the orientation of the staggered magnetization and the 
ground state energy is equal to 



E" 



Considering small deviations from this state, by intro- 
ducing 



together with the corresponding time derivative 
ni w (n;,i,0,ri;,2) , 

we findSi 

L 



(7) 



H.2 



1=1 



-nf 2 (Fc + 7^) - 7^ {ni^ini+i^i + ni^2ni+i,2) } , 
which after Fourier transformation |(2Jl becomes 

-n^in,,! [D^ +B^+-f^l + cosg)] 

-"-g,2»^g,2 [Fc + 7^ (1 + COS q)] } . 

With the time dependence Q , the equations of motion 
are decoupled and we easily deduce the spin-wave spec- 
trum (restoring the in-plane momentum as explained at 
the beginning of Sec. 



ujj = c^k^ + + B^ +-f^ {1 + cosq) , 
u;2 ^c^/c^ + Fc + 7^(1 + cos g) . 



(8) 
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FIG. 1: (Color online). Field- dependence o] the spin-wave 
gaps for different directions of the uniform magnetic field, la- 
beled by the corresponding axis. Ai can be identified as the 
in-plane gap, while A2 corresponds to the out-of-plane gap at 
zero field. These plots are based on the parameters derived in 
Sec.\T]\ I.e., D = 2.5 meV , y/Tl = 5 meV , 7^ = b{meVf, 
= 2.1, and g±_ = 2.5. For a field B > -B^op applied along 
the c axis, we show the curves for momenta q = and q = n. 
A comparison with recent experiments is presented in Fig. 



In the absence of a magnetic field, we find the usual 
in- and out-of-plane gaps at momentum (k, g) — (0, tt), 
see Eq. For an increasing field applied along the a 
axis, the in-plane gap hardens, while the out-of-plane gap 
remains constant. The field-dependence of the spin-wave 
gaps is presented in Fig. ^ 



B. Magnetic field along b 

For a magnetic field applied along the b axis, i.e., B 
Bei,, the Lagrangian reads 



1=1 

[{nff + (nff'j + 2BDni - ^I'^ni ■ n^+i j . (9) 



In this case, the magnetic field first leads to a rotation 
of the n-field around the a axis, as shown in Fig. |21 The 
angle is given by (see also Ref. 0) 



BD 



smc 



+ 272 - 52 ' 
and the ground state energy is equal to 



(10) 



^ 2 



7 



[BDf 



Tc + 272 - 52 



At -Bcrit = ^ ~ 2 IT, it becomes energetically favorable 
to have the largest component of the staggered magne- 
tization directed along the a axis, see Eq. 0. For even 
larger fields, the order parameter is then further rotated 
towards the c axis, with an angle 9 between n and the a 
axis equal to (see also Ref. 0) 



sin 6 



BD 



In this case, the ground state energy is given by 

\2 



(11) 



7^ + - 



[BDY 



Tc + 272 - D2 



Since n2 = 1, we find that for a field of around B k, 90T, 
the staggered magnetization is completely aligned along 
the c direction. 



1. B < B^^t 

Let us now turn to the field-dependence of the spin- 
wave spectrum. In the rotated coordinate system, we 
can express the components of n on even and odd sites 
as 



nil 



ni = 



l-n2 



nf^ cos I 



V 



1 



nf , sm( 



' — 71;, 2 sm ( 
n;,2 cos 9 



''1,1 ~ "-1,2 
Its time-derivative is thus given by 

/ n'l^i \ 

ki^ (-l)'+^n/,2sin6l 



Using the above decomposition of the staggered mag- 
netization, together with the expression for the angle of 
rotation IjlUII . we find the LagrangianSi 



C 



E 

9 



-nl^nq.i [D''-B''+-f^{\ + cosq)\ 



-n*q 2nq,2 [(Fc - +7^ (1 + cosq)) cos 26* + BDsin^ 
-|-2i?cos6' {n* ihg^2 - nq,2nq-i)} ■ 
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FIG. 2: A magnetic field along the b direction leads to a ro- 
tation of the staggered magnetization first around the a axis 
and then around the b axis. 



Note that in contrast to the previous case fSec. lIIlXj) . the 
two modes are now coupled through the last term in the 
above equation. Substituting the time dependence l@J, 
the Euler-Lagrange equations read 



—Lo'^+S^l 2iujB COS. 
-2iLuB cos 9 — Lu 




= 



where 



= [r^ - + 7^(1 + cos g)] cos 261 + sine . 

From here, one finds the spin- wave spectrum for magnetic 
fields B <D 



^(e? + e2 + 4s'cos2 



±V(^i + cos2( 



(12) 



The minimum of the spin-wave gaps is found at (k, q) — 
(0, tt) and their dependence on the magnetic field is plot- 
ted in Fig. U] 



its time derivative being equal to 

Hi = -ri;_i 

\ ni^2 COS0 J 

After Fourier transformation Q and substitution of the 
angle 6 ((TT|l . the Lagrangian © readsSi 



(Fc - I?^ + 7^(1 + COS q)) COS 2^ + 51? sin 6i | 



For fields B > D, the spin- wave dispersions are therefore 
given by 



ujI = c^k^ + - _^ 7^ (1 + cos(?) , 

ul = c^k^ + [Fc - +7^ (1 + cosq)] cos 26* 



BD sine 



(13) 



In the presence of moderate magnetic fields, such that 
cos 20 > 0, which roughly corresponds to B < 65 T, the 
minima of the spin- wave gaps are found at (k, q) ~ (0, tt). 
As soon as 6 > 7r/4, i.e., for extremely high fields, a 
change in momentum of the out-of-plane mode from tt to 
is favorable. 

As illustrated in Fig. ^ the in-plane gap softens for 
B < D and eventually vanishes exactly at the critical 
field. By further increasing the magnetic field, the gap 
hardens again, quickly becoming larger than the zero- 
field value D. In contrast, the out-of-plane gap first 
hardens and then jumps at the transition to decrease 
very slowly for fields B > D. It eventually vanishes in 
the fully polarized state. 



C. Magnetic field along c 

Lastly, we study the effects of a uniform magnetic field 
parallel to the c axis, i.e., B = Bcc- The Lagrangian 
then takes the form.— 



2. -B > -Bcrii 



For magnetic fields B > D, the staggered magnetiza- 
tion is contained in ac plane. Following the same ap- 
proach as below the critical field, we expand the order 
parameter as 



(<-'>' (A 



ni 



"■1,2' 

-nil 



ni 2 sm ( 



Ha 



smt 



nL2 cos( 



1=1 



+B' 



2BDn'' -2 



ni- m+i 



} ■ (14) 



We would like to point out that this case has been ana- 
lyzed earlier in Ref.|2. Nevertheless, we reconsider it here 
to facilitate comparison with experiments in next section. 
For a magnetic field applied along the c axis, there is 
competition between the interlayer coupling term, which 
leads to an antiparallel alignment of the staggered mag- 
netization (still oriented along the h axis) and the term 



6 



linear in the magnetic field, for which a parallel orienta- 
tion is energetically favorable. The mechanism of the ex- 
perimentally observed spin-flop transition follows there- 
fore naturally from the above Lagrangian. 

Let us first calculate the critical field at which the spin- 
flop transition takes place. By comparing the energy of 
the parallel 



El 



2 



(7^ 



with the antiparallel configuration 



El = -^{—l^ + B^ + 2BD) , 



we find 



5flop = « 14 T , 

using the value of 7^ given by Eq. This estimation 
is pretty close to the most recent experimental determi- 
nation of the critical field of around 11.5 T, see Ref. [Til 
but substantially larger than the measurement by Gozar 
et alM, (6.5 T). 

Below the spin-flop transition (for antiparallel align- 
ment), we use the expansion of the n- field given 
in Eq. Q. In momentum representation, the La- 
grangian (|14|l then reads^i 



7"* cos^ q 



Fc + + 7^ ± v/S2D2 +74 cos^q | , 



with two branches, because the momenta are limited to 
the reduced Brillouin zone, i.e., q E [— 7r/2, 7r/2). After 
restoration of the in-plane momentum, we find the spin- 
wave spectrum for fields B < B^^^ 



Lol = C^P + + 7^ ± y/B^D^+Y^CO^ , 

c^k^ + rc + B^ + -f^ ± y^B^D^ + 74 cos2 q . (15) 



The minima of the gaps are therefore located in the lower 
branch at momentum (k,o) = (0,0). These gaps have 
been found earlier in Ref. I3. 

Above the spin-flop transition (for parallel alignment 
of the order parameter), we omit the oscillating factor 
in the expansion of the staggered magnetization Q and 
find the Lagrangian^i 



q 

-nl^nq.y [D^ + BD 
-n;_2"<2,2 [Fc + B'^ + 



- 7^ (1 - cc 
BD - 7^ (1 
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FIG. 3: (Color online). Field-dependence of the in-plane 
spin-wave gaps for different directions of the applied magnetic 
field. Symbols represent the experimental observations shown 
in Fig. 2 of Ref .\^ while the lines correspond to our theoretical 
calculations based on the parameters |J7| ) and the gyromag- 
nettc tensor Note that these values are slightly different 
from those used in Fig. 



The excitation spectrum is thus equal to (see also Ref. 0) 

-7^(1- cos g) , 

7^^ (1 — cos (7) . (16) 



2 2 1 2 

uji = c k 



F, 



52 



BD 



Strictly speaking, the minimum of the gaps is found at 
(k, g) = (0, tt). This implies that there is a change in 
momentum of the low-lying excitations at the spin-flop 
transition. However, the comparison with experiments 
presented in Fig. 13 suggests that Raman scattering actu- 
ally follows the excitation with momentum q = Q. 



IV. COMPARISON WITH EXPERIMENTS 

Let us now compare our findings with recent measure- 
ments by Gozar and co-workers^, who reported the field 
dependence of the spin-wave gaps in LCO observed in 
Raman spectroscopy. As explained in detail in Ref. 123 . 
the in-plane mode can be observed for all three direc- 
tions of the field, whereas the out-of-plane gap is only 
detectable when the magnetic field is applied along the 
orthorhombic h direction. At this stage, it seems appro- 
priate to comment on the uncertainties related to the 
basic parameters that enter the expressions for the spin- 
wave gaps, namely the anisotropics (gaps at zero field) 
and the interlayer coupling. The plots presented in Fig.Q] 
are based on the generally accepted values for these pa- 
rameters. However, their precise values vary from one 
experiment to another, presumably due to slightly differ- 
ent samples. From the experiment under consideration^, 
we deduce 



D ^ 2.15 mey 
7 ^ -uflop 



2 - = 2 {meVf 



(17) 
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While the value of the in-plane gap is in very good agree- 
ment with Ref . 1161 the huge difference in 7^ is quite puz- 
zling. Nevertheless, in Fig. |3| we have reproduced the 
measurements of the in-plane gaps presented in Fig. 2 of 
Ref. ^ together with the results of our theoretical calcula- 
tions based on the parameters given above and our esti- 
mates of the gyromagnetic ratios © . The agreement be- 
tween theory and experiments is quite remarkable. From 
the out-of plane mode observed in Fig. 4 of Ref. 0, we 
conclude that a field applied along the b axis leads to a 
slight hardening of the gap. This is also in good agree- 
ment with our result shown in Fig. ^ 

V. CONCLUSION 

To conclude, we analyze the evolution of the spin- wave 
spectrum in La2Cu04 in a uniform magnetic field. We 
find that a field applied along the orthorhombic b direc- 
tion leads to a rotation of the staggered magnetization 
first in the be and then in the ac plane, until the order 
parameter is completely aligned along the c direction. In 
this case, the field-dependence of the spin-wave spectrum 
is given by Eqs. (|12|l and Ijl^fl . In contrast, for a per- 
pendicular field, there is an ordinary spin-flop transition 
induced by the competition between the magnetic field 
and the interlayer coupling, which leads to the spin- wave 
modes given by Eqs. H15|) and (|16|) . In the case where 
the magnetic field is directed along the a axis, the stag- 



gered magnetization remains aligned along the b direction 
and the evolution of the spin-wave dispersion is given by 
Eq. A comparison with recent measurements of the 
field-dependence of the spin-wave gaps shows a beautiful 
agreement between theory and experiments. We hope 
that this work can stimulate new neutron scattering ex- 
periments that could reconcile the very different spin-flop 
transition fields measured in experiments^iii. 

After completion of this work, we discovered two very 
recent preprints by Benfatto et al.^^''^^ which among 
other questions address the same problem as we do in this 
paper. We would especially like to draw the attention to 
Ref. [23. which contains new measurements of the out- 
of-plane gap and an insightful description of the Raman 
response. Comparing our results, we find that they are 
almost identical, except in the case where a strong mag- 
netic field is applied along the b direction, see Eq. l(T^ . 
While Benfatto and Silva Neto^^ predict a constant out- 
of plane gap for B > D, our results shown in Fig. ^ 
indicate a slow softening of the gap, which eventually 
vanishes at the point where the order parameter is com- 
pletely aligned along the c axis. 
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